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Abstract 

In this paper we give a multiresolution construction in Bergman space. The successful 
application of rational orthogonal bases needs a priori knowledge of the poles of the transfer 
function that may cause a drawback of the method. We give a set of poles and using them we 
will generate a multiresolution in A 2 . We give sufficient condition for this set to be sampling 
sequence for the Bergman space. The construction is an analogy with the discrete affine wavelets, 
and in fact is the discretization of the continuous voice transform generated by a representation 
of the Blaschke group over the Bergman space. The constructed discretization scheme gives 
opportunity of practical realization of hyperbolic wavelet representation of signals belonging to 
the Bergman space if we can measure the values of the transfer function on a given set of points 
inside the unit disc. Convergence properties of the hyperbolic wavelet representation will be 
studied. 
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1. Introduction 

The plan of this paper is as follows. First we introduce a discrete subset of the Blaschke 
group. We will give sufficient conditions for this discrete subset to be sampling set for the 
Bergman space. We present some basic results connected to the Bergman space, we give the 
definition of the voice transform generated by a representation of the Blaschke group on A 2 . 
Using the discrete subset of the Blaschke group we construct a multiresolution decomposition in 
A 2 . First the different resolution spaces will be defined using nonorthogonal basis which shows 
the analogy between the discrete hyperbolic wavelets in A 2 and the discrete affine wavelets in 
L 2 (R). Applying the Gram-Schmidt orthogonalization we consider the rational orthogonal basis 
on the n-th multiresolution level V„. This system is the analogue of the Mamquist-Takenaka 
system in the Hardy spaces, possesses similar properties and is connected to the contractive 
zero divisors of a finite set in Bergman space. We prove that the projection operator P n f(z) 
on the resolution level V n is convergent in A 2 norm to /, is interpolation operator on the set 
the {Jl- Ah, where Ak is defined by (2.7) with minimal norm and P„f(z) — > f(z) uniformly 
on every compact subset of the unit disc. 

1.2. The Blaschke group 

Let us denote by 

(1.1) BJz) := 6 Z ~ b (z G C,a= (b,e) G B := D x T,bz 1) 

1 — bz 



the so called Blaschke functions, where 

(1.2) ©+ := © := {z € C : \z\ < 1}, T := {z € C : \z\ = 1}. 

If a G B, then B a is an 1-1 map on T,and D respectively. The restrictions of the Blaschke 
functions on the set D or on T with the operation (B Sll o _B a2 )(z) := B ai (B a2 (z)) form a 
group. In the set of the parameters B :— D x T let us define the operation induced by the 
function composition in the following way: B S11 o B a2 = _B ai0a2 - The group (B, o) will be 
isomorphic with the group ({_B a ,a G B},o). If we use the notations &j := (bj,€j), j G {1,2} 
and a := (b, e) =: ai o a 2 , then 

(1-3) b= = ■B ( _ i , 2 , 1) (6ie 2 ) ! e = ei = fl (-^ w )4 

l + 0i6 2 £2 1 + e 2 t>ib 2 

The neutral element of the group (B, o) is e := (0, 1) G B and the inverse element of a — (b, e) G B 
is a -1 = (— be, e). 

The integral of the function / : B — > C, with respect to the left invariant Haar measure m 
of the group (B, o) can be expressed as 

(1.4) / B /(a)^(a) = ^// Dl ^ dM6 ^ 

where a = (b, e ,;t ) = (6i + ib 2 , e il ) G D x T. 

It can be shown that this integral is invariant with respect to the left translation a — > ao o a 
and under the inverse transformation a — > a -1 , so this group is unimodular. 

2.1 Special discrete subsets in B and their sampling property 

The one parameter subgroups 

(2.1) Bi :={(r,l) :r G (-1,1)}, B 2 := {(0, c) : e G T} 
generate B, i. e. 

(2.2) a=(0,e 2 )o(0,e 1 )o(r,l)o(0,ei) (a=(re lj e 2 ), re [0, 1), e x , e 2 G T). 

Bi is the analogue of the group of dilation, B 2 is the analogue of the group of translation (see 
[25]). 

The group operation (r, 1) = (n, 1) o (r 2 , 1) in Bi can be expressed using the tangent hyper- 
bolic and its inverse (ath) in the following way 

(2.3) r = I1±IL = thfath n + ath r 2 ) (n, r 2 G (-1, 1)). 

1 + rir 2 

Let denote r = tha, rt = thai, * = 1> 2. Then by 

(n, 1) o (r 2 , 1) = (th Ql , 1) o (tha 2 , 1) = (th {a x + a 2 ), 1), 

it follows that (Bi,o) is isomorphic to (R, +). It is known that (Z, +) is a subgroup of (H, +), 
then Bi = {(thfc, 1), k G 1} is an one parameter subgroup of (Bi, o) (see |27|). 
Let a > 1, denote by 

(2.4) B a = ((r fc ,l): r k = a * ~ ^ , fc G Z 

It can be proved that (B3, o) is another subgroup of (B, o), and (r k , 1) ° (r n , 1) = (ru+n, 1). The 
hyperbolic distance of the points r k , r n has the following property: 



(2.5) p(rfc , r „ ):= ^4^ 

1 - r k r n 



1 



a*"— a "> a"- — a 



Let N(a,k),k > 1 an increasing sequence of natural numbers, N(a, 0) := 1, and consider the 
following set of points zoo '■= 0, 

(2.6) A={zM = V i T, ^ = 0,l,...,iV(a,fe)-l, fc = 0,l,2,...} 
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and for a fixed k £ N let the level k be 



(2.7) Ah = {z k i = r k e x — ,£e {0, 1, N(a, k) - 1} }. 

The points of .4 determine a similar, decomposition to the Whitney cube decomposition of the 
unit disc (see for ex. [24] pp.80). There are two differences: 

1. The radius of the concentric circles are connected to the Blaschke group operation, which is 
important from the point of view when we generate a multiresolution: (r k , l)o(r n , 1) = (r k + n , 1). 
This property is analogue with the property of the dilatation when we generate affine wavelet 
multiresolution levels. 

2. The second: we can choose a and N = N(a, k) such that A will be a set of sampling 
in the Bergman space. In the case of the Whitney cube decomposition the concentric circles 
are defined by r% = 1 — i < r k and we dived this in 2 fe equal parts, but the upper and lower 
density of this set it is not so easy to handle. 

First we will study the following questions: for which choice of a and JV = N(a, k) 

1. will be A uniformly discrete, 

2. will be A an e-net set for some < e < 1, 

3. will be A sampling sequence for Bergman spaces A v l 

To answer these questions let start with some basic definitions and results. For detailed 
exposition concerning these results see for example in [8], |26| . [9]. 

Recall that if z = x + iy 6 D then the normalized area measure is dA(z) = — dxdy. For 
< p < oo, an analytic function / : D — > C belongs to the A p if 



(2.8) / \f(z)f dA(z) < oo. 

JV 

For p — 2 the set A 2 is a reproducing kernel Hilbert space, which is called the Bergman space. 
The reproducing kernel of A 2 is given by the formula 

(2-9) K(z. w) = lr -L_. 

The pseudohyperbolic metric is denned by 

y-z 



l-yz 

A sequence of points V = {z k } of points in the unit disc is uniformly discrete if 

5(r) = inf p{z,,z k ) = 5 > 0. 

For < e < 1, a sequence of points F = {z k : k £ N} of points in the unit disc is said to be e-net 
if each point z G D has the property p(z,z k ) < e for some z k in F. An equivalent statement is 
that 

D = lsC=iA(z fc ,e), 

where A(z k ,e) denotes a pseudohyperbolic disc. 

A sequence of points F = {z k : k £ N} of points in the unit disc is sampling sequence for A p , 
where < p < oo, if there exist positive constants A and B such that 

oo 

(2-10) A\\f\\ p <J2\f(^)\ P (l-\^\ 2 ) 2 <B\\f\\ P , ft A*. 

k=l 

For p = 2 this is equivalent with 

oo 

(2.H) All/H 2 <^|(/,^)| 2 <B||/|| 2 , feA 2 , 

k=l 

where <fi k (z) = K(z, z k )/\\K(z, z k )\\ are the normalized Bergman kernels with w = z k in (2.9). 
These kernel functions are not mutually orthogonal, so far no sequence of distinct points z k £ D 
do the normalized kernel functions form an orthonormal basis. However, this last inequality 
shows that these functions will constitute a frame for A 2 if and only if T — {z k : k £ N} is 
a sampling set for A 2 . A main difference between the Hardy space and the Bergman space is 
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that there is no counterpart of sampling sequences in Hardy spaces. The Bergman spaces A p do 
have sampling sequences, but examples are not so easy to construct. Some explicit examples 
are due to Seip, Duren, Schuster, Horowitz, Luecking (see for ex in [H]). For our purpose we 
need a sampling sequence connected to the Blaschke group operation. 

An A p sampling sequence is never an A p zero-set. A total characterization of sampling 
sequences can be given with the uniformly discrete property and upper and lower density of the 
set see [8], but these densities can be quite difficult to compute. Duren, Schuster and Vukotic 
in [5] gave for sampling sufficient conditions based on the pseudohyperbolic metric, that are 
relatively easy to verify, i.e: for < p < oo, if T is a uniformly discrete e-net with 

1 

e < - 



1+ ' 2 



then is a sampling set for A p . 

Schuster and Varolin [26] improved these sufficient condition. They showed that every uni- 
formly discrete e-net sequence with 



(2.12) e< 1 P 



p + 2 



is sampling set for A p . This last sufficient condition will be used to answer our last questions. 

Theorem 2.1 Let a > 1, and (N = N(a, k),k > 1) a sequence of increasing natural numbers 
and consider the set of points 



.- 2-f 



A = {z kl =r k e*— , i = 0,1,..., N(a,k) -1, k = 0, l,2,...oo}, 

defined as before. Suppose that there exists a = lim^-joo N(a, k)a~ 2k . 

1. If (N(a,k)a- 2k ,k > 1) is increasing sequence and a is finite, then A is uniformly discrete 
and the separation constant satisfies 

S > min | n 



VI + a 2 

2. If (N(a, k)a~' 2k , k > 1) is decreasing and < a < oo, then there exists eo G (0,1) for 
which the set A is eo-net . 

3. If N{a, k)a~' 2k — a, is constant for k > 1, 0<a<oo and 

(a-a- 1 ) 2 + 7v 2 ^-<2p, 

CL 

then A is a sampling set for A p . 

Proof 1. We need to consider two types of situations. The pair of points lie on different 
circles, or they may lie on the same circler of radius r k . Suppose first that the points z k i, z mn 
lie on two different circles of radius r k and r m . Then the generalized triangle inequality for the 
pseudohyperbolic metric (see [8] pp. 38) implies that 

P\Zk£, Zmn I m — n \>n> 0. 

1 - r k r m 

Next suppose that the pair of points lie on the same circler of radius r k , and f / n, then the 
least pseudohyperbolic distant is attained when I = n + 1, then 

2tti I I 2 27Ti I ~ 1 

p(zkt,z kn ) > rk 1 - e N 1 - r k e N 



2r k sin [(1 - r k ) 2 + Art sin 2 -^j 



= {l+[(l-r£)/(2r fc sin(7r/iV))] 2 } 
But sin(7r/iV) > (2/tt)(tv/N) = 2/N, so we deduce that 

p(z ke ,z kn ) > {l + [(l-r 2 )iV/(4r fe )] 2 } 



-1/2 
-1/2 



-1/2 
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We observe that 



(1 - r 2 k )N/(4r k ) = 2fc - _ 2 = jVa- 2fc [l/(l - a~ ik )l 

and p(zki,Zk n ) has a positive lower bound if a — limk^ao N (a, k)a~ 2k < oo and p(z k i,z kn ) > 
— 7==- Combining the two lower bounds we obtain the stated result for the separation constant. 

2. For given z — re* 6 £ D take k and j £ {0, 1, • • • N(a, k) — 1} such that r k < r < r k +i, 



(1 - rr k ) 2 + 4rr k sin 2 6 —^- _ (r - r k ) 2 + 4rr k sin 2 

1 i 7~-\ o \ /1 o\ — 



l-^( 2 , 2fcj ) (l-H)(l-rf) (l-r2)(l-r^) 

(r-r fc ) 2 + 4rr fc ^ _ (q-q^f (q 2fc + 2 - q- 2 ^ 2 )(q 

1 + (l-H)(l-r 2 ) - 1 + 4 + 4 TV 2 " - K(a ' ^ 

If (N(a, k)a~ 2k , k > 1) is decreasing and a = lim^oo JV(a, k)a~ 2k € (0, oo), then the last term 
is upper bounded by 

t, , {a -a' 1 ) 2 a 2 2 
K:=l+ i , — i- + -^tt 2 . 

4 4a 2 

Then for 

eo = %/l - T/K, 

we have p(z, z k j) < eo- 

3. Using (2.12) we have that, if (a — a -1 ) 2 + 7I " 2 ^ < 2p, then eo < \J ~fp2 > consequently A is 
sampling sequence for A p . 

Remark From this theorem we obtain that if A is a sampling set for A p , then 

(a — a -1 ) 2 < 2p 

therefore a must be in the interval (1, v^p+V^p +I ) Then we can always choose N = N(a,k) 
big enough, such that the the sampling condition will be satisfied. From the point of view of 
computations and to have on every circle the less possible numbers, for p — 2 a convenient 
choice is a = 2, and N(2,k) = 2 2k+fi for k > 1 with /3 a fixed integer. Then a = 2 fi , and the 
smallest value for j3 for which the sampling condition is satisfied is f3 = 3, then on the fc-th circle 
we will have iVi(2, fc) = 2 2fc+3 equidistant points corresponding to the roots of order 2 2k+i of 
the unity. For a = V2 for sampling we need N\(\f2, k) — 2 k+2 points. 

2.2 The continuous voice transform on Bergman space 

Our goal is to construct a multiresolution analysis based on the set A in the Bergman space. 
First we will need a few results connected to the Bergman space [8] and the voice transform of 
the Bergman space. 

The set A 2 = A 2 (D) with the scalar product 



(2.14) {f,g):= / f(z)g(z) dA(z) 

JO 

is a Hilbert space. An analytic function in the unit disc of the form 

oo 

/(*):=£<**" (, 6 D) 

n=0 

belongs to the set A 2 if and only if the coefficients satisfies 

oo ^ 

^ |C " |2 (^TT) <0 °- 

n=0 V ' 

The Bergman space _A 2 (B) is a closed subspace of L 2 (D). 
For each z£D the point-evaluation map 

r z :A 2 (D)->C, T ,(f) = f(z) 
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is a bounded linear functional on A 2 (D). Each function / G A 2 (3) has the property 

i/(*)i^- l/2 *(*r l ii/iu»(D) (*ed), 

where 5{z) = dist(z,T). From this it follows that the norm convergence in A 2 (B) implies the 
locally uniform convergence on D. Therefore, by the Riesz Representation Theorem there is a 
unique element in A 2 (D), denoted by K{.,z), such that 

/(*) = r,(/) = (f,K(.,z)), (/ G A 2 (B),z G D). 

The function 

K-.HxH^C with K(.,z) £ A 2 (V>) 

is called the Bergman kernel for D. Taking /(£) = K(£,w) for some i«6Dwe conclude that 
the kernel function has the following symmetry property: 



K(z,w) = - / w)K(£, z)d^d$ 2 =K{w,z). 
^ Jo 

For any orthonormal basis {<Pj,j = 0, 1, 2, ...} for j4 2 (B) one has the representation 

oo 

(2.15) K(£,z) = J2<Pi(0w(z), K,«)6Dx», 

3 = 1 

with uniform convergence on compact subsets of D x 1. The set of functions 



{ipi(z) = ^/(j + l)z j , z G D, j = 0, 1, 2, ...} 
form an orthonormal basis in A 2 (H), consequently 

(1 - z£) 2 

The explicit formula for the kernel function shows that 

1 f „^ 1 _ t 2, 



(2.11) /(*) = - / f(t)——-d£ 1 db (/ e A 2 (B), z G D). 

^ Jo (l-£z) 2 

Applying this formula in particular for f(z) = (1 — fz)~ 2 for fixed z in the disc we obtain that 

1 



ii^")ii 2 = ^ / 1 1 f 



U 4 dtid& (1 " M 2 ) 2 



Kfz,«) > 0. 



The voice transform on Bergman space is induced by a unitary representation of the Blaschke 
group on the Bergman space. Results connected the voice transform on Bergman space were 
published in [22] . 

Let consider the following set of functions 



2.10 F a (z) := ^ 1 J 6 ^ (a = (6,e) GB.0GB). 

1 — bz 

This function induces a unitary representation on the space A 2 . Namely let define 

U,f -[F^ffoB- 1 (a G B, / G A 2 ), 
then this is a representation, i.e.: 

i) U x .y = U x oU y (x, y G B), 

ii) B 3 x — }■ U x f G A 2 is continuous for all / G A 2 . 

And more ?7 a (a G B) is an unitary, irreducible square integrable representation of the group 
on the Hilbert space A . 
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The voice transform of / G A 2 generated by the representation U and by the parameter 
g G A 2 is the (complex-valued) function on B defined by 

(Y B f)(x) := (f,U m g) (x G B, /, g G A 2 ). 
In |22] a direct proof of the analogue of Plancherer formula was given, i.e.: 
[V P1 f,V P2 g}=4ir( Pl ,p 2 ) (f,g) (f,g, P i, P 2 G ^ 2 (D)), 

where 

[F, G] := [ F{a)Gja)dm(a), 

m is the Haar measure of the group B. 

This transform is in same relation with the Blaschke group and the Bergman space as the 
affine wavelet transform with the affine group and the L 2 (R) (see 25 ). Indeed let consider 
(G, o) equal to the the affine group, where 

G = {£( a ,b) :i->l:(n,l))6l*X R}, 

t(a,b)( x ) = ax + b, R* := R\ {0}, £1 o £ 2 (x) = £i(£ 2 (x)) = a ia2 x + ai b 2 + bi. 
The representation of the affine group G on L 2 (R) is given by 

U( a ,b)f(x) = H~ 1/2 f{a~ l x - b), 

where a is the dilatation parameter, and b the translation parameter. 

The continuous affine wavelet transform is a voice transform generated by this representation 
of the affine group: 

W*f{a,b) = \a\- 1/2 f f(W(a-H-b)dt={f,U (a , b) il>), f,i/>€ L 2 (R). 
Jm 

There is a rich bibliography of the affine wavelet theory (see for example [4] [6] 1151 1171 1181 1191 
[33]. One important question is the construction of the discrete version, i.e., to find if) so that the 
discrete translates and dilates 

VVfc = 2- n/2 4>(2- n x - k) 

form a (orthonormal) basis in L 2 (R) and generate a multiresolution (see [7 , 9, 17, 19, 20, 21] 
etc.). The construction of the discrete version is connected to the discrete subgroup of the affine 
group, generated by the following set: 

G n , k = {^(a-»,-k) : K -> K : n G Z, k G Z}. 

The discretization of the voice transform can be achieved using the unified approach of the 
atomic decomposition elaborated by Feichtinger and Grochenig 11 . This general description 
can be applied when the integrability condition of the voice transform is satisfied. In a recent 
paper [23] it is shown that, the integrability condition in the Bergman space it is not satisfied. 
This motivates to find discrete multiresolution decomposition in Bergman spaces? 

2.3 Multiresolution analysis in the Bergman space 

We start with the general definition of the affine wavelet multiresolution analysis in L 2 (W). 

Definition 2.2.1. Let Vj, j G Z be a sequence of subspaces of L 2 (R). The collections of 
spaces {Vj, j G Z} is called a multiresolution analysis with scaling function <j> if the following 
conditions hold: 

1. (nested) Vj C Vj+i 

2. (density) iWj = L 2 {R) 

3. (separation) (~]Vj = {0} 

4- (basis) The function cj> belongs to Vo and the set {2 n ^ 2 cf)(2"x — fc), k G Z} is a (orthonor- 
mal) bases in Ki- 
ln multiresolution analysis, one decomposes a function space in several resolution levels and 
the idea is to represent the functions from the function space by a low resolution approximation 
and adding to it the successive details that lift it to resolution levels of increasing detail. 
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Wavelet analysis couples the multiresolution idea with a special choice of bases for the 
different resolution spaces and for the wavelet spaces that represent the difference between 
successive resolution spaces. If V„ are the resolution spaces Vo C Vi C .... C V n —, then the 
wavelet spaces W n are defined by the equality W n © V n = Vn+i- 

In the construction of affine wavelet multiresolutions the dilatation is used to obtain a higher 
level resolution (f(x) G V n <4> f(2x) £ V^+i) and applying the translation we remain on the 
same level of resolution. This field has also a rich bibliography (see for example [Tll2ll6ll4l[7l ll0l 
Q^QTlQBEnilllllO]). Using the subgroup B 3 of the Blaschke group, a discrete subgroup of B 2 
and the representation we give a similar construction of the affine wavelet multiresolution in the 
Bergman space. To show the analogy with the affine wavelet multiresolution we first represent 
the levels V„ by nonorthogonal bases and then we construct an orthonormal bases in V„ and 
give also an orthogonal basis in W n which is orthogonal to V n . We will show that in the case 
of this discretization the analogue of the Malmquist-Takenaka systems for Bergman space, will 
span the resolution spaces and the density property will be fulfilled, i.e., Ufcli ^fe ~ ^ m norm. 
Similar multiresolutin results for the Hardy space were obtained by the author in |21] , 

We show that the projection P n f on the n-th resolution level is an interpolation operator in 
the unit disc until the n-th level, which converges in A 2 norm to /. 

Let consider a > f, denote by r k ~ , k G N, TV = N(a,k) a sequence of natural 

numbers such that a — N(a, k)a~ 2k satisfies 

2 

n ^ I -1\2 . 2 O, 

< a < oo, (a ~ a ) + 7r — < 4. 

or 

Let us consider the set of points 

A = {z u =r k e iS 7T, £ = 0,1,...,N-1, k = 0, 1, 2, ...}, 
and for a fixed k G N let the level k be 

(2.17) A h = {z M = r k e l2 Tr, £ G {0, 1, AT — 1} }. 

Due to Theorem 2.1 A is a sampling set for A 2 . This implies that the set of normalized kernels 
(ip k i(z) — jY^==^yj, k = 0, 1, ■ • • , £ — 0, 1, ■ ■ ■ TV — 1) will constitute a frame system for A 2 . 
From the frame theory [12] or from atomic decomposition results (see Theorem 3 of [37]) follows 
that every function / from A 2 can be represented 

f( z ) = c kl (p kl (z) 

(k,l) 

for some {c k t } G £ 2 , with the series converging in A 2 norm. The determination of the coefficients 
it is related to the construction of the inverse frame operator (see Q2]), which is not an easy task 
in general. This is the reason why we try to construct other approximation process for / G A 2 
such that the determination of the coefficients follow an exactly defined algorithmic scheme. 

Let us consider the function </?oo = 1 and let Vo = {c, c G C}. 

Let us consider the nonorthogonal hyperbolic wavelets at the first level 

(2.18) y M (s) = (U (zUil) -tp )(z) = (1 5 , £ = 0,l,--.JV(o,l). 

v " ' (1 — zuzy 

They can be obtained from 931,0 using the analogue of translation operator which in the unit disc 
is a multiplication by a unimodular complex number, and from 950,0 using first the representation 
operator U( Tli i)-i followed by the translation operator: 

(2.19) V\A Z ) = <Pi,o(ze~ 1 ^ r > ) ) = (^(r 1 ,i)-i^o,o)(2e _ ^T ) ). 

Let us define the first resolution level as follows 
(2.20) 

JV(a,l)-l 

Vi = {/ : D — > C, f(z) = c ,o<A>,o + J2 c ^i,i, coo, c M £ C, I = 0, 1, • • • TV(a, 1) - 1.}. 

£=0 
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Let us consider the nonorthogonal wavelets on the n-th level 

(2.21) <p n , t {z) = (U (Zne<1) -ip )( Z ) = (1 ~_^ ] , e = 0,l,-,N(a,n)-l, 

\ i Z n gZ) 

which can be obtained from <p n ,o using the translation operator, and from 930,0 using the repre- 
sentation U(t r _ 1 ,i)o(ri i)) _I > ano - the translations 

(2-22) Vnit (z) = ([/( {r „_ 1 ,i)o(, 1 ,i))-iPo)(ze" lw ^ ) ). 

Let us define the n-th resolution level by 

n N{a,k)-1 

(2.23) V n = {/ : D -> C, f(z) = £ C M</?M> c k ,e G C }. 

fc=o e=o 

The closed subset V n is spanned by 

(2.24) {<p hi t, £ = 0,1, N(a, k) - 1, fc = 0, n}. 

Continuing this procedure we obtain a sequence of closed, nested subspaces of A 2 for z£D 

(2.25) V C Vi C V 2 C K C .... 

Due to Theorem 2.1 the normalized kernels {ipti{z) = t^==^ , k = 0, 1, • ■ ■ , i = 0, 1, ■ ■ • N(a, fc) — 
1} form a frame system for A 2 which implies that this is a complete and a closed set in norm, 
i.e., 

(2.26) |J V n = A 2 , 

new 

consequently the density property it is satisfied. 

For a = 2 and N{2,n) = 2 2 "+ 3 , if a function / G V„, then (7 (ri ,i ) -i/ G K+i. For this it is 
enough to show that 

P(n,i)-'WW = ^(n.D-iK^.D-iPo)]^"'^)) = 

(2.27) [(^(r, +1 ,i)-iPo)](ze" < 2 2 ^+«+ 3 ) G V n+1 , k = 1, .., n, £ = 1, ...2 2fe+3 - 1. 

From now on for simplicity we will deal with this case, but in the general case we can always 
choose N(a, k) such that the previous condition to be valid. 

Since the set A is a sampling set it follows that is a set of uniqueness for A 2 , which means 
that every function / G A 2 is uniquely determined by the values {f(zke)}- In the paper of Kehe 
Zhu [38] described in general how can be recaptured a function from a Hilbert space when the 
values of the function on a set of uniqueness are known and developed in details this process in 
the Hardy space. At the beginning we will follow the steps of the recapturation process but we 
will combine this with the multiresolutin analysis. The set 

(2.28) ( * £ = 0,l,...,2 2fc+3 -l, fc = 0,l,...,n.} 
{{1- z kl zy J 

is a nonorthogonal basis in V n . 

Using Gram-Schmidt orthogonalization process they can be orthogonalized. Denote by ijjk,t 
the resulting functions. They can be seen as the analogue of the Malmquist -Takenaka system 
in the Hardy space. This functions can be obtained using the following two methods. The first 
arises from the orthogonalization procedure. To describe this let reindex the points of the set A 
as follows 01 = zoo, 0,2 = Z10, a-i = 211,04 = 212, 033 = 21,31,034 = z 2 ,o, ■ ■ • a m = zu ■ • ■ , k = 
0, 1, I = 0, 1, 2 2fc+3 - 1, and denote by K(z, z ht ) = n J_ := K(z, a m ) 

<f>oo(z) = <p(ai, z) 
</>ke(z) = 0(01,02, 
a m = Zki, m > 2. 
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K(z,ni) 

|I#(-,<h)|| ' 

,o m , z) = K(z,a„ 



Em — 1 
i=i 1 



(Oi,02, 



_ ,\ (K(.,q m ),(/>(ai,a2,...,n i ,.)> 

■y u " z ) ||V>(ai,a 2 ,...,a i ,.)|| 2 : 



(1) 



Thus the normalized functions {ipkt(z) = ^■jrjp & = 1, 2, ••• ,£ — 0, 1, ■ ■ ■ 2 2fc+3 } became an 
orthonormal system. Applying similar construction in Hardy space we get in this way the 
Malmquist-Takenaka system. They can be written in a nice closed form using the Blaschke 
products. Unfortunately in our situation this is not the case and the properties of the system 
can be not seen from the previous construction. 

Another approach is given by Zhu in [5S] . He shows that the result of the Gram-Schmidt 
process are connected to some reproducing kernels and the contractive divisors on Bergman 
spaces. Let denote A m = {ai, ai, ■ ■ ■ a m }. Let H Am the subspace of A 2 consisting of all functions 
in A 2 which vanish on A m . H Am is a closed subspace of A 2 and denote by K Am the reproducing 
kernel of H Am ■ These reproducing kernels satisfies the following recursion formula 

ts t \ v i \ K Am (z,a m+1 )K Am (a m+1 ,w) 
K A m+1 (z,w) = K Am (z,w) r ,m > 0, 

K-A m (dm+l, a m + i ) 

(2.29) K Ao := K(z,w) = - 1 

(1 — WZJ 

The result of the Gram-Schmidt process can be expressed as 

K(z,ai) K Al (z,a 2 ) K Am _ 1 (z,a m ) 



y / K(ai,a 1 )' \J K Al {a2,a 2 y y / K Am _ 1 (a m , a m ) ' 

Then 

(2.30) i>kt(z) = == 

\/AA m _ 1 (a m , dm) 

and is the solution of the following problem 

sup{ J Re/(a m ) : / G H Am _ t , ||/|| < 1}. 

This extremal functions in the context of the Bergman spaces have been studied extensively in 
recent years by Hedenmalm [13]. The main result in [13] is that the function 

K Am _ ± (z, am) 



\/K Am _ 1 (a m ,a m ) 

is a contractive divisor on the Bergman space, vanishes on A m -i, and if A is not a zero set 
for A 2 , as is in our case, the functions converge to as m -> oo. In Hardy space the partial 
products of a Blaschke product corresponding to a nonzero set own all these nice properties. 
From the Gram-Schmidt orthogonalization process it follows that 

(2.31) V n = span{i; k , e , £ = 0,1, 2 2k+z - 1, k = 0~^}. 

The wavelet space W„ is the orthogonal complement of V n in V n +i- We will prove that 

(2.32) W n = span{i> n +i,t, 1 = 0, 1, 2 2n+5 - 1}. 
Indeed, every function / £ A 2 can be recovered using the, 

(2.33) /W = j f_Z^_ <M(w) . 

If / £ V n , one has f(z) = X]fe=o X^=o + _1 c k,e i Pk,i C A 2 , then using (2.30) we obtain that 

n 2 2fc + 3 -l 

(V'n+lj,/) = ^2 c k,i{<l>n+l,j,<Pk,l} = 

fc=o e=o 

n 2 2fc + 3 -l 

Ck ^ - rl)^ n+1> e(z h i) = 0, j = 0, 1, ...2 2n+5 - 1. 

fc=0 f=0 

We have proved that for f £ V n 

(2.34) (f,ip n+1 , j )=0, 
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which is equivalent with 

(2.35) Vn+u-LVn, (j = 0,l,...,2 2 " +5 -l). 
From 

(2.36) V n +i = Vn@span{tp n +i, j , j = 0,l,...,2 2n+5 -l} 
it follows that W„ is an 2 2 ( n+1 ' +3 dimensional space and 

(2.37) W n = span{ip n+ i,t, £ = 0, 1, 2 2n+5 - 1}. 

Summary. We have generated a multiresolution in A 2 and we have constructed a rational 
orthogonal wavelet system which generates the levels of the multiresolution. 

2.4 The projection operator corresponding to the n-th resolution level 

Let us consider the orthogonal projection operator of an arbitrary function / G A 2 on the 
subspace V n given by 

n 2 2fc + 3 -l 

(2.38) P n f(z) = J2 E </>^mWm(*)- 

fc=0 t=o 

This operator is called the projection of / at scale or resolution level n. 

Theorem 2.2 For f € A the projection operator P n f is an interpolation operator on the 
points zm — rke l 2 2k + 3 , (I — 0, 2 2fc+3 — 1, k = 0, n), is norm convergent in A 2 to f, i.e. 

\\f - Pnf\\ ->0, n->oo, 

uniformly convergent inside the unit disc on every compact subset, and is the solution of minimal 
norm interpolation problem. 

Proof Let consider N = 1 + 2 5 + • • • + 2 2n+3 and the corresponding kernel function of the 
projection operator 

n 2 2fc + 3 -l 

k n (z,o = J2 E mz)m*) = 

fc=0 t=0 



(2 39) V K A m -i(z,a m ) ( K Am _ 1 {j, a n l ) \ = K Am _ 1 (z, a m )K Am _ 1 (am, Q 
From the recursion relation (2.29) it follows that 

JV 

(2.40) k n (z, o = E o - Ka ^ (*■ 0) = ^(*> e) - e) 

m— 1 

From this relation it follows that the values of the kernel- function in the points zte, (I 
0, 2 2fc+3 - 1, k = 0, n) are equal to 

(2.41) K(zu,Z) ' 



(l-*«0 2 

Every function f £ A 2 can be recovered using the Bergman projection 
Therefore 

(2.42) P n f( Zke )=[ - !M - dA(w) = f(z mj ) 0' = 0,....,2 2m+3 -l, m = 0,...,n). 
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We obtain that P n f is interpolation operator for every / G A 2 on the set Um = o-4m- 

Because of 2.26 and 2.31 {ipk,t, ft = 0,oo, £ = 0, 1, 2 2fe+3 - 1} is a closed set in the 
Hilbert space A 2 , we have that that ||/ — P n f\\ — > as n — > oo. Since convergence in A 2 norm 
implies uniform convergence on every compact subset inside the unit disc, we conclude that 
P™f( z ) ^ fi z ) uniformly on every compact subset of the unit disc. From Theorem 5.3.1 of [24] 
there exists a unique /„ g V n with minimal norm such that 

(2.43) f n (z mj ) = f(z mj ), 0=0,....,2 2m+3 -l, m = 0,...,n), 

/„ is uniquely determined by the interpolation conditions and is equal to the orthogonal pro- 
jection of / on Vn, thus f n {z) = P n f{z). 

2.5 Reconstruction algorithm 

In what follows we propose a computational scheme for the best approximant in the wavelet 
base {^ kA , £ = 0,1, -, 2 2fc+3 - 1, ft = 0, n}. 

The projection of / £ A 2 onto V n +i can be written in the following way: 

(2.44) P n+1 f = P n f + Q n f, 
where 

2 2 "+ 5 -l 

(2.45) Qnf(z):= J2 (f,i>n+l,e)i>n+lAz)- 

1=0 

This operator has the following properties 

(2.46) Qnf(z ke ) = 0, k = 1, n, < = 0, 1, 2 2n+3 - 1. 

Consequently P n contains information on low resolution, i.e., until the level An, and Q n is the 
high resolution part. After n steps 

n 

(2.47) Pn+if = Pif + Y,Q"f- 

k=i 

Thus 

(2.48) VWi = ^o0Wo0Wi0...0W„. 

The set of coefficients of the best approximant P n f 

(2.49) {b ke = (f,iP Ke ), £ = 0.1,...2 2k+3 ~1 fc = 0,l,...,n} 

is the (discrete) hyperbolic wavelet transform of the function / £ A 2 . Thus it is important to 
have an efficient algorithm for the computation of the coefficients. 

The coefficients of the projection operator P n f can be computed if we know the values of 
the functions on {Jl^ Ak- For this reason we express first the function tpk.e using the bases 
i}Pk' ,i' V = 0, 1, ...2 2fc +3 — 1, k' = 0, ft), i.e. we write the partial fraction decomposition of 
tfrkl ■ 

k-l 2 2k '+ 3 -l 1*1 

(2.50) *w = £ E c k ,, t , +Y,*°j {1 ^ r 

k'=o i'=o K s; 2=0 v JS; 

Using the orthogonality of the functions (ipk',i' £' = 0, 1, ...2 2fc +3 — 1, ft' = 0, ft) and the 
reconstruction formula 

fc_l 2 2k '+ 3 -l i 

(2.51) SknS(m = {lpnm,1pki) = E E C k' ,£' i>n.m(Zk> e>) + ^ C k,j1pn,m(Zkj) , 

k'=0 t'=0 j=0 

(m = 0, 1, ...2 2n+3 - 1, n = 0, ft). 
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If we order these equalities so that we write first the relations (2.51) for n = k and m = 
t,i-\, ...,0 respectively, then for n = k - 1 and m = 2 2(fc - 1)+3 - 1, 2 2(fc - 1)+3 - 2, ...,0, etc., 
this is equivalent to 



(2.52) 



o 





W 



( ipkA z k,e) 

4>k,i-2(Zk,e) i>k,l-2(Zk,l-l 

V ipoo(zk,e) 4>oo{zk,e-i) ipoo{z h ,i-2) 



\ 








/ \ 




Ck,e-i 




Ck,e-2 




\ coo / 



This system has a unique solution (cs~I, Ck,e-i, Ck,e-2, •••, c oo) T - If we determine this vector, 
then we can compute the exact value of (/, tpk,e) knowing the values of / on the set Ufe=o *^k- 

Indeed, using again the partial fraction decomposition of tpk,e and the reconstruction formula 
formula we get that 

k — 1 2 2fc —1 £ 



fc-1 2 2k ' -1 I 

E ^k r 7'}{zk',l')+^Ck~jf{Zk,j)- 
3=0 



(2.53) 



k'=0 £'=0 



Summary. Measuring the values of the function / in the points of the set A = Ufc=o -^ fe ^ ® 
we can write the operator (P n f, n £ N) which is convergent in A 2 norm to /, is the minimal 
norm interpolation operator on the set the Ufc=o anc ^ P™f( z ) f( z ) uniformly on every 
compact subset of the unit disc. 
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